Disappearance of static magnetic order and evolution of spin fluctuations in Fei+^Se^^Tei- 
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We report neutron scattering studies on static magnetic orders and spin excitations in the Fe-based chalco- 
genide system Fei+^Sej^Tei-i, with different Fe and Se compositions. Short-range static magnetic order with 
the "bicolhnear" spin configuration is found in all non-superconducting samples, with strong low-energy mag- 
netic excitations near the (0.5, 0) in-plane wave-vector (using the two-Fe unit cell) for Se doping up to 45%. 
When the static order disappears and bulk superconductivity emerges, the spectral weight of the magnetic ex- 
citations shifts to the region of reciprocal space near the in-plane wave-vector (0.5, 0.5), corresponding to the 
"collinear" spin configuration. Our results suggest that spin fluctuations associated with the collinear mag- 
netic structure appear to be universal in all Fe-based superconductors, and there is a strong correlation between 
superconductivity and the character of the magnetic order/fluctuations in this system. 

PACS numbers: 74.70.Xa, 75.25.-j, 75.30.Fv, 61.05.fg 



I. INTRODUCTION 

Since the discovery of the first high temperature super- 
conductor in the 1980's, there has been a continuing ef- 
fort to understand the origin of high-Tc superconductivity. 
Studies on the cuprate systems seem to suggest that there is 
an intimate relationship between superconductivity and mag- 
netism,'^ and recently this has been shown to be the case 
also for the newly discovered Fe-based superconductor fam- 
ilies.^ '" With non-superconducting parent compounds that 
have static antiferromagnetic (AF) order, charge doping in the 
Fe pnictides gradually suppress the AF order and induces su- 
perconductivity, for both the RFeAsO ("1111" system)""'* 
and AFe2As2 ("122" system) families i'^i'^ In addition, a "spin 
resonance" has been observed in the "122" system by inelastic 
neutron scattering, '^"-^ showing a sharp increase of the mag- 
netic scattering intensity at the "resonance" energy when the 
system goes into the superconducting phase. For the pnic- 
tides, both the static magnetic order in the parent compound 
and the "resonance" in the superconducting compounds occur 
around the in-plane wave vector (0.5,0.5) (using the 2-Fe unit 

cell), suggesting a "collinear" or "C-type" AF structure [see 
Fig.[I](a)]^12i25 

The situation is slightly different in another Fe-based su- 
perconductor family, the iron chalcogenide Fei+^SexTei-^r 
(the "1:1" compound). Here the parent compound has a "bi- 
collinear" or "E-type" AF order [see Fig. [I](b)],^^^' modu- 
lated along the (0.5,0) in-plane direction. Furthermore, the 
substitution of Se for Te, which induces superconductivity, 
does not directly modify the density of electrons in the con- 
duction bands. Despite these differences, the spin resonance 
observed in the superconducting compositions of this family 
has been found to occur at the same (0.5, 0.5) in-plane wave- 
vector^^— as in the Fe pnictides, but rotated 45° from the 
ordering wave-vector of the parent FeTe compound. This sug- 



gests that the superconducting mechanism in the Fe pnictides 
and chalcogenides are likely to be quite similar Nevertheless, 
there has also been a report implying that superconductiv- 
ity can coexist with the "bicoUinear" structure on the atomic 
scale. ^^ Therefore, understanding how the magnetic structure 
(static or dynamic) evolves from the E-type bicollinear config- 
uration of the parent compound to the C-type collinear config- 
uration in the superconducting region in the 1 : 1 system is an 
important problem. There have been a number of theoretical 
studies^'''^* that provide some insight on this issue, but clearly 
a more systematic experimental study is highly desirable. 

In this paper, we present our work using neutron scattering 
to probe the magnetic order/fluctuations in a few samples from 
the 1 : 1 family for Se dopings ranging from 30% to 50% and 
with varying superconducting properties. Our results suggest 
that static magnetic order exists in all non-superconducting 
samples. (Here, by non-superconducting we mean an ab- 
sence of bulk superconductivity.) This order is short-ranged 
and occur at in-plane wave-vectors of the type (0.5, 0). For 
the fully superconducting samples, no static magnetic order 
is found. With the disappearance of static magnetic order, 
the associated low energy magnetic excitations near (0.5, 0) 
also go away, as one might expect. Magnetic excitations 
near (0.5, 0.5) gradually become dominant as the material be- 
comes more superconducting. While Se doping plays an es- 
sential role, it is clearly not the only determining factor re- 
garding the superconductivity and the magnetic correlations. 
Samples with similar Se doping but differing in Fe content 
can have very different superconducting properties and corre- 
sponding magnetic structures/fluctuations. Our results clearly 
indicate that static bicollinear (E-type) magnetic orders in the 
1 : 1 system compete with and suppress superconductivity. Su- 
perconductivity only appears when the system evolves toward 
a fluctuating collinear C-type magnetic correlations, which 
appear to be universal across all known Fe-based supercon- 



TABLE I: List of the Fei+^Sei^Tei-i samples used in our measure- 
ments, withi tiieir composition (S, x), superconducting transition tem- 
perature (Tc), room-temperature lattice parameters (from powder x- 
ray diffraction), and crystal mass. 



Sample 


S 


X 


T, 


a 


c 


mass 








(K) 


(A) 


(A) 


(g) 


SC30 


0.00 


0.3 


14 


3.815 


6.140 


12.7 


NSC30 


0.05 


0.3 


- 


3.808 


6.120 


7.4 


SC50 


0.00 


0.5 


14 


3.811 


6.129 


9.0 


NSC45 


0.05 


0.45 


- 


3.807 


6.047 


6.4 



at around 14 K, SC50 is clearly better in quality as far as su- 
perconducting volume fraction is concerned. With a consid- 
erable portion of its bulk volume being non-superconducting, 
it is possible that there is phase separation in SC30. In fact, 
when measuring different small pieces {^ 1 mm size) from 
the same SC30 sample, the superconducting volume can vary 
from < 10% to ^ 80%, suggesting that the superconduct- 
ing and nonsuperconducting phases could be macroscopically 
separated in this sample. The other samples, NC30 and NC45, 
are mostly non-superconducting, with no more than 1 % of the 
volume giving a superconducting response. 



ductor families. 



II. EXPERIMENT 
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FIG. 1 : (Color online) (a) Schematic of the "collinear" C-type AF 
spin structure, with scattering intensities mainly around (0.5,0.5). 
The square shows a unit cell with two Fe atoms, (b) Schematic of 
the "bicollinear" E-type AF spin structure, which contribute mainly 
to scattering intensities near (0.5, 0). (c) The schematic diagram of 
the neutron scattering measurements in the (HKO) zone. Dashed 
lines denote linear scans performed across (0.5, 0.5) and (0.5, 0) in 
the text. 



The single-crystal samples used in the experiment were 
grown by a unidirectional solidification method. The sam- 
ples, their nominal compositions, and various characteristic 
properties are listed in Table I. The bulk susceptibility results 
in Fig. |2] were obtained using a superconducting quantum in- 
terference device (SQUID) magnetometer. From the magne- 
tization measurements we can see that, although both super- 
conducting samples show evidence of diamagnetic response 
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FIG. 2: (Color online) ZFC magnetization measurements by SQUID 
for SC30 (red), NSC30 (blue), SC50 (black) and NSC45 (green). The 
inset shows the same data from the non-superconducting samples 
with different scale. 



Neutron scattering experiments have been caiTied out on 
the triple-axis spectrometers SPINS (inelastic scattering mea- 
surements of magnetic excitations for SC30 and NSC30, and 
all elastic measurements for static magnetic order), and BT- 
7 (inelastic scattering measurements for SC50 and NSC45) 
located at the NIST Center for Neutron Research. We used 
horizontal beam collimations of Guide-80'-S-80'-240' (S rep- 
resents "sample") for the inelastic scattering measurements 
on SPINS with fixed final energy of 5 meV and a cooled Be 
filter after the sample to reduce higher-order neutrons; col- 
limations of Guide-open-S-80'-240' were used for the elas- 
tic measurements on SPINS, with an additional Be filter be- 
fore the sample. At BT7, we used beam collimations of 
open-50'-S-50'-240' with fixed final energy of 14.7 meV and 
two pyrolytic graphite filters after the sample. The inelastic 
scattering measurements have been performed in the (HKO) 
scattering plane, as indicated in Fig. [itc). The data are de- 
scribed in reciprocal lattice units (rl.u.) of (a*,6*,c*) = 
(27r/a, 27r/6, 27r/c). The elastic scattering measurements 



have been taken in the (HOL) scattering plane instead, since 
the order in the parent compound occurs at half integer L 
values (AF order along the c-axis). All data have been nor- 
malized into absolute units (/is^-eV^^/Fe), using incoherent 
elastic scattering intensities from the samples. 



in. RESULTS AND DISCUSSION 
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FIG. 3: (Color online) Elastic neutron scattering measurements per- 
formed on SC30 (top) and NC45 (bottom) near (0.5, 0, 0.5). (a) and 
(b) are intensity profiles along [100] direction (//-scans); (c) and (d) 
are scans along [001] direction (L-scans). (e) and (f) show the mag- 
netic peak intensity at (0.5, 0, 0.5) vs. temperature. Corresponding 
scans measured at T = 34 K are used as background, and have 
been subtracted from all the data shown. The error bars represent the 
square root of the number of counts. 



The static long-range magnetic order in the parent com- 
pound FeTe appears near (0.5, 0, 0.5), suggesting a bicollinear 
E-type AF structure. With small Se doping, it was suggested 
that the order should gradually become short-ranged^^'"*^ and 
eventually disappear. Our results, however, suggest that the 
order can still remain with relatively large Se dopings. 

We performed elastic magnetic scattering measurements on 
all samples. For the SC50 sample, there is no elastic magnetic 
intensity at (0.5, 0, 0.5), while magnetic peaks are observed 
for all three other samples. In Fig. [3] we plot H and L scans 
through this AF wave-vector for the SC30 and NSC45 sam- 
ples at r = 4 K. The same measurements for NSC30 have 
been published in a previous report (see Fig. 2 of Ref. uTD . 
The H and L scans performed at higher temperature (T=34 K) 
show no peak structure and are therefore used as backgrounds 
to be subtracted from the data. All peaks are much broader 
than the resolution, indicating the short-ranged nature of the 



magnetic order The if -scans are peaked near but not exactly 
at H = 0.5, similar to the results reported for NSC30 and an- 
other sample with 27% Se doping.'^ The L-scans, however, 
are qualitatively different. In previous reports on lower Se 
doped 1:1 compounds, -^'''^ the L-scan peaks around L — 0.5, 
and intensity always goes to zero at L = 0. So there the mag- 
netic order is always AF along the L-direction, whether short- 
or long-ranged. Here we see that after background subtrac- 
tion, the scattering intensity at L = is still appreciable. This 
suggests that although the magnetic order still has a modu- 
lation along the L-direction, which peaks around L = 0.5, 
favoring an AF configuration between Fe planes, the order 
has become much more two-dimensional. In the NSC45 sam- 
ple, the 3D long-range bicollinear AF magnetic order of the 
parent compound has not been entirely destroyed, but rather 
greatly reduced to 2D short-range order The ordered moment 
per Fe is 0.015(4) /i^, much less than the value in the parent 
compound with long-range order. ^^ It is nevertheless, enough 
to destroy superconductivity. With this static magnetic order 
present, even with 45% Se doping, bulk superconductivity is 
still not achieved. In the SC30 sample, although the sample 
does show a superconducting phase transition at around 14 K, 
the superconducting volume is smaller than for the SC50 sam- 
ple. The ordered moment is about 0.006(2) /i^/Fe, also much 
less than that in the NSC45 sample, indicating that this order 
may be coming from only part of the sample. Therefore there 
is likely a phase separation in this SC30 sample, where two 
phases, one superconducting and another one with a short- 
range magnetic order, coexist. 

With the tendency of forming static bicollinear magnetic 
structures in the non-superconducting samples, it is natural to 
expect to see magnetic excitation spectra around the (0.5,0) 
in-plane wave-vector as well. Previous work has shown that 
the energy dispersion and intensity modulation along the L 
direction for magnetic excitations in the 1:1 compound is 
small. ^^'^^ We can therefore choose to perform the inelastic 
scattering measurements in the {HKO) plane for L = 0. In 
Fig.m we plot our results taken near (0.5, 0, 0). The top two 
panels [Fig. |4] (a) and (b)] show energy scans at (0.5, 0, 0) 
at T = 4 K and 25 K. Measurements for NSC45 and SC50 
were taken on BT7 with a relatively coarse energy resolution 
(FWHM - l.VmeV) compared to those on SPINS (NSC30 
and SC30, FWHM ~ 0.3 meV), and have a large, resolution- 
limited tail from scattering at hu ~ 0. Constant-energy scans 
atHuj = 0.5, 2, and 5 meV [Fig.|4](c) to (h)] along the K di- 
rection across (0.5,0,0) clearly show that for NSC30, SC30 
and NSC45, there is significant spectral weight at low energies 
here. For both 30%-Se samples, where we have measurements 
with higher energy resolution, one can see that the intensity at 
huj — 0.5 meV increases on warming from 4 K to 25 K. The 
increase is much less pronounced at fiw = 5 meV. This be- 
havior is likely due to a transfer of spectral weight from the 
elastic peak into low energy channels when the static order 
dissolves with heating. For SC50, the low-energy spin exci- 
tations near (0.5, 0, 0) are weak, and not strongly temperature 
dependent, which is consistent with the fact that there is no 
static order near (0.5,0) in this sample. The two small peaks 
near K — ±0.5 observed from samples SC30 and SC50 at 
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FIG. 4: (Color online) Magnetic excitations for Fei+^Se^iTei-a; 
measured around (0.5, 0,0). The left and right columns show the 
magnetic peak profiles for lowest temperature and 25 K respectively, 
(a) and (b) Constant-Q scans at (0.5, 0, 0) taken at low T and 25 K. 
(c-h) Constant-energy scans at (0.5, K, 0) at Ikv =0.5, 2 and 5 meV. 
A fitted if-independent background has been subtracted from all data 
sets. The error bars represent the square root of the number of counts. 
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FIG. 5: (Color online) Magnetic excitations for Fei+^Se^^Tei-j, 
measured around (0.5, 0.5, 0). The left and right columns show the 
magnetic peak profiles for lowest temperature and 25 K respectively, 
(a) and (b) Constant-Q scans at (0.5, 0.5, 0) taken at low T and 25 K. 
(c-h) Constant-energy scans at (0.5,0.5,0), taken along the trans- 
verse direction [as shown in Fig.[T](a)] at fiu = 5, 6.5 and 12 meV. 
A fitted constant background has been subtracted from all data sets. 
The error bars represent the square root of the number of counts. 



Hljj = 5 meV, suggest that there is additional spectral weight 
developing near the (0.5,0.5) wave-vector, corresponding to 
dynamic collinear spin correlations in the superconducting 
samples. 

In Fig. |5] we show measurements near (0.5,0.5,0). For 
SC50, a clear "resonance" is observed when comparing the 
energy scans performed at 4 K and 25 K [see Fig. |5](a) and 
(b)]. In panels (c) to (h), constant-energy scans at hu =5, 
6.5 and 12 meV performed in the direction transverse to Q = 
(0.5,0.5) are shown. Similar to SC50, but less pronounced, 
we can also see a "resonance" feature in SC30 in the scans of 
hu! = 6.5 and 5 meV. 

For the non-superconducting samples, there is no tempera- 
ture effect observed for data taken between 4 K and 25 K. For 
NSC30, we did a constant-energy scan near (0.5, 0.5) only at 
hu! — 6.5 meV, and the intensity is very low compared to ei- 
ther its own magnetic scattering near (0.5, 0) or those from 
the other samples near (0.5, 0.5). It is clear that the low en- 
ergy spin excitations are mostly focused around (0.5,0) for 
NSC30. The NSC45 sample has Se doping very close to 
SC50, and also very similar magnetic excitation spectrum near 
(0.5, 0.5) compared to that from the latter in its normal state 
(T — 25 K). However, with no superconducting transition, its 
spectrum at low temperature (T = 4 K) does not differ much 
from that at T = 25 K. 

The implications of our results are very clear for NSC30 
and SC50. For NSC30, a short-ranged static magnetic order 



is present at low temperature near (0.5, 0), corresponding to 
a 3D bicollinear E-type spin structure. Its low energy mag- 
netic excitations are also focused near (0.5, 0). With the static 
order present, no superconductivity is achieved in this sam- 
ple. For the SC50 sample, there is no static order and the low 
energy magnetic excitation spectrum is mostly shifted to the 
(0.5, 0.5) region, corresponding to collinear C-type spin cor- 
relations. Similar to the situation in the "122"'^''^ or "1111" 
systems,"""' this collinear configuration without static order 
appears to favor superconductivity. 

The results for NSC45 are more complicated. Here, with 
Se doping close to SC50, the magnetic excitations near the 
(0.5, 0.5) point are rather similar to the superconductor, ex- 
cept that the "resonance" feature is missing. Apparently, hav- 
ing magnetic excitations near (0.5, 0.5) associated with the 
collinear spin configuration is not sufficient for superconduc- 
tivity to emerge. Also having a static 2D-like magnetic order 
with bicollinear structure is able to completely suppress su- 
perconductivity in this sample. Of course, the tetragonal crys- 
tal structure gives no energetic distinction between the order- 
ing wave vectors (0.5, 0) and (0.5, 0.5), so that the magnetic 
configuration is relatively soft. There is likely a mixed phase 
where the bicollinear and collinear magnetic configurations 
coexist, on a microscopic level, similar to that of the mixed 
C-E phase in manganites,^^ as suggested in Ref.l24l 

The case for SC30 is, in fact, quite intriguing. A 2D-like 



short-range static order exists at low temperature, while low 
energy magnetic excitations are found both around (0.5, 0) 
and (0.5,0.5) with comparable spectral weight. Therefore, 
the magnetic excitation spectrum actually looks very simi- 
lar to that in NSC45, yet there is bulk superconductivity in 
SC30 when the static magnetic order is also present. Com- 
pared to NSC45, the ratio of spectral weight near (0.5,0.5) to 
that near (0.5,0) is higher in SC30, indicating a larger vol- 
ume of the sample favoring a coUinear spin configuration. 
The "resonance" occurs below T^, showing an enhancement 
of spectral weight only near (0.5,0.5). This indicates that 
superconductivity only exists in the part of the sample with 
dynamic collinear spin correlations. Although it is conceiv- 
able that the static order and superconductivity could coexist 
in the same domains as suggested by previous /^jSR work,-'^ 
it is also possible to have a system with macroscopic phase 
separation, where the volume of local collinear or bicollinear 
regions are large enough to form separate domains. In this 
case, the features near (0.5, 0) (elastic magnetic peak and low 
energy magnetic excitations), and those near (0.5, 0.5) come 
from different regions. This scenario would be consistent with 
the (varying) susceptibility results for different small pieces 
taken at different locations from this sample, and agrees with 
results from all other samples where static magnetic order and 
superconductivity do not coexist locally. 

Why would samples with similar Se content (e.g., NSC30 
vs. SC30, NSC45 vs. SC50) show dramatically different be- 
haviors? It is clear from previous work that the Fe content 
has a significant impact on the magnetism.^^-^^ For example, 
higher Fe content in the parent compound can drive the or- 
der from commensurate to incommensurate,^* while its effect 
has been less clear for the superconducting region. It seems 
unlikely that the excess Fe atoms act simply as isolated mag- 
netic moments that destroy the superconductivity, since our 
observations of variations in static magnetic order and low 
energy spin excitations cannot be explained in such a simple 
manner. Theoretically, it has been predicted that lowering the 
height of the chalcogen (Te/Se) positions can drive the "1:1" 
system from the bicollinear to the collinear spin configura- 
tion,^'*-'^ and Fe interstitials will certainly have an impact on 
Fe-Te/Se-Fe bond lengths and bond angles. There are several 



reports concerning the effect of excess Fe on the lattice param- 
eters'*"^^; our results indicate that the lattice parameters a and 
c both decrease slightly with increased Fe content (holding 
Se constant), which does not appear to be entirely consistent 
with the trends reported by others. To make further progress 
on this issue, it will likely be necessary to characterize the 
microstructure associated with specific compositions.^^Sdl 



IV. SUMMARY 

Despite uncertainties in the underlying cause of differences 
among our samples, it is evident that the magnetic struc- 
ture/fluctuations are intimately correlated with the supercon- 
ductivity in the 1 : 1 system. If the magnetic correlations in the 
system favor a bicollinear (E-type) spin configuration, which 
sometimes can eventually lead to static order, superconductiv- 
ity is suppressed. The collinear (C-type) spin configuration, 
which is universal across the superconducting regions of all 
Fe-based superconductor families, is necessary, but not suffi- 
cient, for the emergence of bulk superconductivity. There are 
cases where magnetic correlations favoring these two config- 
urations coexist and compete. Overall, our results suggest that 
magnetic correlations are important in the Fe-based supercon- 
ductor families, and the proper tuning of these correlations 
may be the key for enhancing superconductivity. 



Acknowledgments 

We thank Weiguo Yin and Wei Ku for useful discussions. 
Work at Brookhaven is supported by the Office of Basic En- 
ergy Sciences, Division of Materials Science and Engineering, 
U.S. Department of Energy (DOE), under Contract No. DE- 
AC02-98CH10886. JSW andZJX ai-e supported by the Center 
for Emergent Superconductivity, an Energy Frontier Research 
Center funded by the U.S. DOE, Office of Basic Energy Sci- 
ences. The SPINS spectrometer at the NCNR is supported 
in part by the National Science Foundation under Agreement 
No. DMR-0454672. 



* J. M. Tranquada, H. Woo, T. G. Perring, H. Goka, G. D. Gu, 
G. Xu, M. Fujita, and K. Yamada, Nature 429, 534 (2004). 

^ H. F. Fong, P. Bourges, Y Sidis, L. P. Regnault, A. Ivanov, G. D. 
Gu, N. Koshizuka, and B. Keimer, Nature 398, 588 (1999). 

^ G. Xu, G. D. Gu, M. Hucker, B. Fauque, T. G. Perring, L. P Reg- 
nault, and J. M. Tranquada, Nature Phys. 5, 642 (2009). 

" H. A. Mook, P Dai, S. M. Hayden, G. Aeppli, T. G. Perring, and 
F. Dogan, Nature 395, 580 (1998). 

^ Y Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. Am. 
Chem. Soc. 130, 3296 (2008). 

' H. Takahashi, K. Igawa, K. Arii, Y Kamihara, M. Hirano, and 
H. Hosono, Nature 453, 376 (2008). 

' X. H. Chen, T. Wu, G. Wu, R. H. Liu, H. Chen, and D. F Fang, 
Nature 453, 761(2008). 



Z.-A. Ren, W. Lu, J. Yang, W. Yi, X.-L. Shen, C. Zheng, G.-C. 

Che, X.-L. Dong, L.-L. Sun, F. Zhou, et al.. Chin. Phys. Lett. 25, 

2215 (2008). 

F-C. Hsu, J.-Y Luo, K.-W. Yeh, T.-K. Chen, T.-W. Huang, P M. 

Wu, Y.-C. Lee, Y.-L. Huang, Y.-Y Chu, D.-C. Yan, et al., Proc. 

Natl. Acad. Sci. USA 105, 14262 (2008). 

K.-W. Yeh, T.-W. Huang, Y.-l. Huang, T.-K. Chen, F-C. Hsu, 

P M. Wu, Y.-C. Lee, Y.-Y. Chu, C.-L. Chen, J.-Y. Luo, et al.. 

Euro. Phys. Lett. 84, 37002 (2008). 

C. de la Cruz, Q. Huang, J. W. Lynn, J. Li, W. R. li, J. L. Zarestky, 

H. A. Mook, G. F Chen, J. L. Luo, N. L. Wang, et al.. Nature 453, 

899 (2008). 

J. Zhao, Q. Huang, C. de la Cruz, S. Li, J. W. Lynn, Y Chen, 

M. A. Green, G. F Chen, G. Li, Z. Li, et al.. Nature Mater. 7, 953 



(2008). 

Q. Huang, J. Zhao, J. W. Lynn, G. F. Chen, J. L. Luo, N. L. Wang, 
and P. Dai, Phys. Rev. B 78, 054529 (2008). 
C.-H. Lee, A. lyo, H. Eisaki, H. Kito, M. T. Fernandez-Diaz, 
T. Ito, K. Kihou, H. Matsuhata, M. Braden, and K. Yamada, J. 
Phys. Soc. Jpn. 77, 083704 (2008). 

J. Zhao, Q. Huang, C. de la Cruz, J. W. Lynn, M. D. Lumsden, 
Z. A. Ren, J. Yang, X. Shen, X. Dong, Z. Zhao, et al., Phys. Rev. 
B 78, 132504 (2008). 

Y Qiu, M. Kofu, W. Bao, S. H. Lee, Q. Huang, T. Yildirim, 
J. R. D. Copley, J. W. Lynn, T. Wu, G. Wu, et al., Phys. Rev. 
B 78, 052508 (2008). 

J. Zhao, D. T. Adroja, D.-X. Yao, R. Bewley, S. Li, X. F. Wang, 
G. Wu, X. H. Chen, J. Hu, and P Dai, Nature Phys. 5, 555 (2009). 
H. Chen, Y. Ren, Y Qiu, W. Bao, R. H. Liu, G. Wu, T. Wu, Y L. 
Xie, X. F Wang, Q. Huang, et al. Euro. Phys. Lett. 85, 17006 
(2009). 

A. D. Christianson, E. A. Goremychkin, R. Osborn, 
S. Rosenkranz, M. D. Lumsden, C. D. Malliakas, I. S. Todorov, 
H. Claus, D. Y Chung, M. G. Kanatzidis, et al.. Nature 456, 930 
(2008). 

M. D. Lumsden, A. D. Christianson, D. Parshall, M. B. Stone, 
S. E. Nagler, G. J. MacDougall, H. A. Mook, K. Lokshin, 
T. Egami, D. L. Abernathy, et al., Phys. Rev. Lett. 102, 107005 
(2009). 

S. Chi, A. Schneidewind, J. Zhao, L. W. Harriger, L. Li, Y Luo, 
G. Cao, Z. A. Xu, M. Loewenhaupt, J. Hu, et al., Phys. Rev. Lett. 
102, 107006 (2009). 

S. Li, Y Chen, S. Chang, J. W. Lynn, L. Li, Y Luo, G. Cao, Z. Xu, 
and P Dai, Phys. Rev. B 79, 174527 (2009). 
D. S. Inosov, J. T. Park, P Bourges, D. L. Sun, Y Sidis, A. Schnei- 
dewind, K. Hradil, D. Haug, C. T. Lin, B. Keimer, et al.. Nature 
Phys. 6, 178 (2010). 

W.-G. Yin, C.-C. Lee, and W. Ku, arXiv:1003.0512vl (2010). 
S. Li, C. de la Cruz, Q. Huang, Y Chen, J. W. Lynn, J. Hu, Y.-L. 
Huang, F-C. Hsu, K.-W. Yeh, M.-K. Wu, et al., Phys. Rev. B 79, 
054503 (2009). 

W. Bao, Y Qiu, Q. Huang, M. A. Green, P Zajdel, M. R. Fitzsim- 
mons, M. Zhernenkov, S. Chang, M. Fang, B. Qian, et al., Phys. 
Rev. Lett. 102, 247001 (2009). 
J. Wen, G. Xu, Z. Xu, Z. W. Lin, Q. Li, W. Ratcliff, G. Gu, and 



J. M. Tranquada, Phys. Rev. B 80, 104506 (2009). 

Y Qiu, W. Bao, Y Zhao, C. Broholm, V. Stanev, Z. Tesanovic, 

Y C. Gasparovic, S. Chang, J. Hu, B. Qian, et al., Phys. Rev. Lett. 
103, 067008 (2009). 

H. A. Mook, M. Lumsden, A. Christianson, B. C. Sales, R. Jin, 
M. A. McGuire, A. Sefat, D. Mandrus, S. Nagler, T. Egami, et al., 
arXiv:0904.2178vl (2009). 

J. Wen, G. Xu, Z. Xu, Z. W. Lin, Q. Li, Y Chen, S. Chi, G. Gu, 
and J. M. Tranquada, Phys. Rev. B 81, 100513(R) (2010). 

D. N. Argyriou, A. Hiess, A. Akbari, I. Eremin, M. Kor- 
shunov, J. Hu, B. Qian, Z. Mao, Y Qiu, C. Broholm, et al., 
arXiv:0911.4713vl (2009). 

S.-H. Lee, G. Xu, W. Ku, J. S. Wen, C. C. Lee, N. Katayama, Z. J. 
Xu, S. Ji, Z. W. Lin, G. D. Gu, et al., arXiv:0912.3205vl (2009). 
S. Li, C. Zhang, M. Wang, H.-q. Luo, E. Faulhaber, A. Schnei- 
dewind, J. Hu, T. Xiang, and P Dai, arXiv:1001.1505vl (2010). 
W. Bao, A. T. Savici, G. E. Granroth, C. Broholm, K. Habicht, 

Y Qiu, J. Hu, T. Liu, and Z. Mao, arXiv:1002.1617vl (2010). 

R. Khasanov, M. Bendele, A. Amato, P. Babkevich, A. T. 

Boothroyd, A. Cervellino, K. Conder, S. N. Gvasaliya, H. Keller, 

H. H. Klauss, et al., Phys. Rev. B 80, 140511(R) (2009). 

C. Fang, B. A. Bernevig, and J. Hu, Euro. Phys. Lett. 86, 67005 

(2009). 

N. Katayama, S. Ji, D. Louca, S.-H. Lee, M. Fujita, T. J. Sato, 

J. S. Wen, Z. J. Xu, G. D. Gu, G. Xu, et al., arXiv:1003.4525vl 

(2010). 

T. Hotta, M. Moraghebi, A. Feiguin, A. Moreo, S. Yunoki, and 

E. Dagotto, Phys. Rev. Lett. 90, 247203 (2003). 

C.-Y Moon and H. J. Choi, Phys. Rev. Lett. 104, 057003 (2010). 

T. M. McQueen, A. J. Williams, P W. Stephens, J. Tao, Y Zhu, 

V. Ksenofontov, F. Casper, C. Felser, and R. J. Cava, Phys. Rev. 

Lett. 103, 057002 (2009). 

T. M. McQueen, Q. Huang, V. Ksenofontov, C. Felser, Q. Xu, 

H. Zandbergen, Y S. Hor, J. Allred, A. J. Williams, D. Qu, et al., 

Phys. Rev B 79, 014522 (2009). 

R. Viennoisa, E.Gianninia, D. v. d. Marela, and R. Cerny, 

arXiv:0911.2081v2(2009). 

R. Hu, E. S. Bozin, J. B. Warren, and C. Petrovic, Phys. Rev. B 

80,214514(2009). 



